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Abstract 

We study time dependence of exchange symmetry properties of Bell states when two qubits 
interact with local baths having identical parameters. In case of classical noise, we consider a 
decoherence Hamiltonian which is invariant under swapping the first and second qubits. We find 
that as the system evolves in time, two of the three symmetric Bell states preserve their qubit 
exchange symmetry with unit probability, whereas the symmetry of the remaining state survives 
with a maximum probability of 0.5 at the asymptotic limit. Next, we examine the exchange 
symmetry properties of the same states under local, quantum mechanical noise which is modeled 
by two identical spin baths. Results turn out to be very similar to the classical case. We identify 
decoherence as the main mechanism leading to breaking of qubit exchange symmetry. 

PACS numbers: 03.65.Yz 03.65. Ta 
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I. INTRODUCTION 



Since the early days of quantum mechanics, it has been known that certain quantum states 
have a mysterious non-local behavior [1]. The phenomenon responsible for these non-local 
correlations among the subsystems of a composite quantum system is called entanglement 
[2]. Quantum entanglement, having no classical counterpart, is believed to be one of the 
characteristic features of quantum mechanics. Besides its foundational importance for the 
quantum theory, entanglement is also considered as the resource of quantum computation, 
quantum cryptography and quantum information processing [3]. In recent years, it has 
been extensively studied with various motivations [4]. However, entanglement of quantum 
systems, as all other quantum traits, is very fragile when they are exposed to external 
disturbances, which is inevitably the case in real world situations. 

Decoherence, the process through which quantum states lose their phase relations irre- 
versibly due to interactions with the environment, is crucial for understanding the emer- 
gence of classical behavior in quantum systems [5]. It also presents a major challenge for 
the realization of quantum information processing protocols since protection of non-local 
correlations against undesirable external disturbances is essential for the reliability of such 
protocols. Consequently, understanding the decoherence effect of the environment on entan- 
gled systems is an important issue. This problem has been currently addressed in literature, 
considering both local and collective interactions of qubits and qutrits with the environ- 
ment. While some authors examined the effects of classical stochastic noise fields [6-9], 
others studied the same problem for large spin environments [10-15]. 

In this work, we focus on a different aspect of a decoherence process of entangled states. 
Certain two-qubit entangled states have the property that they remain unchanged under 
the exchange of two qubits. We will concentrate on a decoherence model which also has an 
exchange symmetry, i.e., having a Hamiltonian invariant upon swapping the first and second 
qubits. Our goal is to understand how the exchange symmetry properties of symmetric pure 
states alter as the quantum system evolves in time for a symmetric Hamiltonian which 
embodies the effect of local and identical noise fields on qubits. More specifically, we will 
investigate the exchange symmetry properties of three of the four Bell states. Bell states 
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are defined as maximally entangled quantum states of two-qubit systems and given as 
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We will only consider the first three of these states which are symmetric under exchange 

operation. However, our discussion can be extended to include anti-symmetric states like 
|i?4). The first three Bell states are among the symmetric pure two-qubit states which can 
be represented in the most general case by the density matrix 
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where |ap + 2|cp + |6p = 1. After classical noise calculations, we will briefly discuss the 
exchange symmetry properties of the same states for local and quantum mechanical noise 
which is modeled via two identical large spin environments. 



II. LOCAL CLASSICAL NOISE 



We assume that the two qubits are interacting with separate baths locally and the initial 
two-qubit system is not entangled with the local baths. The model Hamiltonian we consider 
was first introduced and studied by Yu and Eberly [6] and can be thought as the represen- 
tative of the class of interactions which generate a pure dephasing process that is defined 
as 

H{t) = ~tJi[nA{t){a, ® /) + nB{t){I ® a,)], (6) 
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where we take h — 1 and cr^ is the Pauh matrix 

/ 1 n \ 

(7) 

Here fi is the gyromagnetic ratio and nA{t),nB{f) are stochastic noise fields that lead to 
statistically independent Markov processes satisfying 

{n,{t)) = 0, (8) 



{niit)n,it'))^^Jit-t'), 
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where (■■■) stands for ensemble average and rj(i=A,B) are the damping rates associated 
with the stochastic fields and nB{t). 

The time evolution of the system's density matrix can be obtained as 

Pit) = {u{t)p{o)uHt)), (10) 

where ensemble averages are evaluated over the two noise fields nA{t) and nB{t) and the 
time evolution operator, U{t), is given by 

U{t) = exp[-i [ dt'Hit') ]. (11) 

JQ 

The resulting density matrix in the product basis {|00), |01), |10), |11)} can be written as 
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P211B P22 P23lAlB P2A1A 

PillA p32lAlB P33 P^AIB 
\ PAllAlB Pi2lA PiZlB PU 

where pij stands for the elements of the initial density matrix, p(0), and ja, 1b are given by 
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For our purposes, we want our two local baths to be identical in a sense that they have the 
same dephasing rate T. Therefore, we let = Fb = F. The resulting density matrix of the 
system with the consideration of identical baths is now given by 
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where 7a = 7a = 7- 



III. OPERATOR-SUM REPRESENTATION OF DECOHERENCE 



To examine the symmetry properties, we need to express the dynamical evolution of p{t) 
in terms of quantum operations. The decoherence process of our quantum system can be 
regarded as a completely positive linear map that takes an initial state p(0) and maps 
it to some final state p{t) [3]. For every completely positive linear map there exists an 
operator-sum representation which is known as Kraus representation [16-18]. The effect of 
the map is given by 
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p(t) = $(p(0)) = 5^ir,(i)p(0)4(^), 
where are the Kraus operators which satisfy the unit trace condition 
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The Kraus operator approach provides an elegant way to study the decoherence process. 
In order to describe the internal decoherence dynamics of the system, all we need to know is 
the Kraus operator set which inherently contains the entire information about environment. 
The operator-sum representation of our completely positive linear map, $(p), which reflects 
the effect of the stochastic process, can be obtained by studying the mapping called Choi- 
Jamiolkowski isomorphism [17-18]. In our investigation, it turns out that the effect of the 
mapping, $(p), on the two-qubit system can be expressed by a set of four Kraus operators 
as 
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where a;(t), a{t) and (3{t) are given by 



(20) 



uj{t) = Vl - oc{t) = 7(t) - 1, I3{t) = 7(t) + 1. 



(21) 



Since different environmental interactions may result in the same dynamics on the system, 
the operator-sum representation of a quantum process is not unique. The collective action of 
our set of the four Kraus operators {Ki, K2, K^, K4} on the density matrix of the two-qubit 

quantum system are equivalent to the collective action of another set of Kraus operators 
{El, E2, E^, E4} if and only if there exists complex numbers Uij such that Ei = ^jUijKj 
where Uij are the elements of a 4 x 4 unitary matrix [3]. The unitary freedom will provide 
us an easy way to introduce exchange symmetry condition. 



IV. EXCHANGE SYMMETRY OF BELL STATES UNDER DECOHERENCE 

In order to analyze the exchange symmetries of symmetric Bell states I-B2) and 
I-B3), we exploit the unitary freedom on the operator-sum representation. Consider the 
most general 4x4 unitary matrix with complex elements 

/ ?y,ii ?y,io ?y,iQ ?y,i,i \ 



Uii U12 M13 Mi4 

U21 U22 U23 U24 

U31 U32 M33 «34 

\ 1*41 1*42 1*43 1*44 / 
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where WU = I. 
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The mapping described by the Kraus operators {Ki, K2, K^, K4} , via unitary freedom, 
is equivalent to the mappings described by the following four Kraus operators 



P T^. r ^^/'l , '^^-"3 , '^«m2 «^m3 , /3^iM4 

^M = Diag[--^ + ^ + ^,--^-^ + ^. 

^ 2 + 2 ' V2 + 2 + 2 

where = 1, 2, 3, 4. 



(23) 



A. Exchange Symmetries of \Bi) and 1^2) 



Having calculated all possible Kraus operator sets, we are in a position to evaluate the 
possible final states when the initial state is \Bi) or I-B2). The density matrices of possible 
final states are obtained as 

E,{t)p^^{0)El{t) E,{t)p^^{0)El{t) 



Tr{E,it)pB^{0)El{t)) 



Tr{E,{t)pB^{0)El{t)y 
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where p = 1,2,3,4 and p^*{0) = \Bi){Bi\ for i = 1,2. The explicit forms of the density 
matrices p^^ {t) and p^^ {t) are given by 
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Obviously, the symmetry condition given in Eq. (5) brings no restriction on these density 
matrices. Thus, it is guaranteed that the Bell states \B-i) and IS2) always preserve their 
exchange symmetry as they evolve in time under our model Hamiltonian. 



B. Exchange Symmetry of IB3) 



The density matrices of the possible final states for {Bs) are written as 

where 11 — 1, 2, 3, 4 and p^^ (0) = \Bs) {Bs]. The exphcit form of the density matrix p^^ (t) is 
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As can be seen from the form of the density matrix of the most general two-qubit sym- 
metric pure state in Eq. (5), for possible final states to be symmetric we need all non-zero 
elements of the matrix in Eq. (29) to be equal to each other, that is, r = s . This condition 
can only be satisfied in case of u^2 = 0. We can immediately conclude that it is impossible 
for all of the possible final states to be symmetric since any 4x4 unitary matrix has to 
satisfy the condition that |wi2p + 1^22^ + |ii32p + |'if42p = 1- Thus, IS3) cannot evolve in 
time under our model Hamiltonian in a way that preserves its qubit exchange symmetry 
with unit probability. In other words, the exchange symmetry of this two-qubit state has 
to be broken with some non-zero probability. Considering the symmetry of the initial state 
and the Hamiltonian this is a very interesting result. A natural question is the maximum 
probability of finding a symmetric possible final state as the system evolves in time. In 
order to answer this question, we need to consider three different cases, namely, the cases of 
having one, two or three symmetric possible final states. 

If we assume only one of the possible final states to be symmetric, say the outcome of Ei 
{ui2 = 0), then the probability of getting a symmetric output state is given by 



Psym(t ^ 00) = - \uis + Uul"^ . 



(31) 



If we assume two of the possible final states to be symmetric, say the outcomes of Ei and 
E2 (mi2 = 0, U22 = 0), then the probability of having a symmetric output state is given by 

Psym{t ^ 00) = ^\ui3 + Uu\'^ + ||m23 + lt24|^- (32) 



Finally, if three of the possible final states are symmetric, say the outcomes of Ei, E2 and 
E3 {ui2 — 0, 1^22 = 0, 1*32 = 0), then the probability of having a symmetric output state is 
given by 

Psymit ^ 00) = ^luis + Uul"^ + ^11*23 + ^^24^ + ^1^33 + U34\'^ . (33) 

In all of these possible cases, the maximum probability of finding a symmetric final state 
turns out to be 0.5. 



V. LOCAL QUANTUM NOISE 

When it comes to modeling the baths as large spin environments, one of the simplest 
decoherencc models, introduced in [19], is that of two central spins interacting with N 
independent spins through the Hamiltonian [13] 

Ni N2 
-H" = Ci2 ^ hiOikCTikz + C-2Z ^ flt02k(^2kz- (34) 
k=l k=l 

This model describes two central spins, with ^-component operators Ci^ and C2z, coupled 
to bath spins represented by ankz, where n — 1,2 labels the baths and k — 1, 2, 3, 
labels the individual spins. All spins are assumed to be 1/2 and cu, c^z and Onkz denote the 
corresponding Pauli matrices. If wc assume that the central spins are not entangled with the 
spin baths at t = 0, the initial state wiU be in product form |^'(0)) = |^'c(0))|^'ai(0))|^' 0-2(0)) 
where 

l^c(O)) = (ottI IT) + a^l Ti) + a^tl it) + «ul (35) 

with 

l*an(0)) = (g)(anfc| \nk) + ^nk\ ink)). (36) 
k=\ 

where | \nk) and | ink) are eigenstates of Onkz with eigenvalues +1 and -1, respectively, and 

The reduced density matrix of two central spins at later times will be given by tracing 
out the bath degrees of freedom from the total density matrix of the system, p(i), as Pc(^) = 
Tr^jpii) where subscript a means that trace is evaluated by summing over all possible nk 
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states and p{t) — \^{t)){^{t)\. The resulting reduced density matrix in product basis 
{|TT),|n),UT),|U)}isfound to be 



Pc 



|2 



2 



f I«ttP 

where the decoherence factors ri(t) and r2(t) are given by 

k=l 



\ 



(37) 



(38) 



In general both expansion coefficients ank, Pnk and interaction strengths cUnk a-re random. 
For our purposes, we will assume that the baths are identical, which means we let expansion 
coefficients and interaction strengths of the two baths be equal to each other as aik = a2k = 
Oik, Pik — P2k — Pk and ujik — uj2k — '^fe- This assumption implies that the decoherence 
factors of two baths are equal so that ri{t) — r2{t) — r{t). Thus, the reduced density matrix 
of two central spins is simplified to 



|2 



|2 1^ 12 



where 



N 

k=l 



(39) 



(40) 



We immediately observe that the form of Pc under the assumption of identical baths is very 
similar to the form of the output density matrix we obtained for classical noise Hamiltonian. 
In particular, when the initial expansion coefficients ak and Pk are equal to each other, we will 
have exactly the same form of the mapping obtained in Sec. II. Hence, decay of r(t) to zero at 
later times and the form of the possible Kraus operators in this case guarantee that the qubit 
exchange symmetry properties of symmetric Bell states \Bi), |i?2) and \B^) interacting with 
two local large spin environments will be the same as their behavior under local stochastic 
noise fields. Since we interpret decay of r(t) as a signature of decoherence, we identify 
decoherence as the main source of spontaneous breaking of qubit exchange symmetry. 
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VI. CONCLUSION 



We examined the time evolution of exchange-symmetric Bell states for local noise Hamil- 
tonians having the same symmetry. For both classical and quantum noise, we found a rather 
unexpected result that not all Bell states preserved their symmetry. In fact, we observed that 
exchange invariance property survived with a maximum probability of 0.5 at the asymptotic 
limit. We conclude that breaking of exchange symmetry for some possible final states is a 
characteristic feature of decoherence. 
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